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Abstract. We provide realistic mock-catalogs of cosmic rays above 40 EeV, 
for a pure proton composition, assuming their sources are a random subset of 
ordinary galaxies in a simulated, volume-limited survey, for various choices of 
source density: 10~ 3 5 Mpc -3 , 10~ 4 Mpc -3 and 10~ 4 ' 5 Mpc -3 . The spectrum 
at the source is taken to be E~ 2 3 and the effects of cosmological redshifting as 
well as photo-pion and e+ e — energy losses are included. 

In order to develop and test analysis tools for ultrahigh energy cosmic ray 
(UHECR) correlation studies it is necessary to have mock catalogs of UHECRs under 
various realistic scenarios for the sources and magnetic deflections. There is no need 
at present to postulate a "new physics" origin for UHECRs, so the sources of UHECRs 
will be associated with the distribution of matter in the universe. The most popular 
types of source candidates reside in ordinary galaxies - for instance, AGNs and GRBs - 
and to first approximation these galaxies are a random subset of all galaxies. Here we 
provide mock UHECR catalogs assuming the sources are randomly selected galaxies. 
In this case, the key characteristic of a source population is its number density. The 
small number of doublets in the Auger events above 57 EeV leads to the lower bound 
of 61 on the number of sources producing them pQ; this in turn leads to a lower bound 
on the source density of about 3 x 10 -5 Mpc" 3 Simulations of the spectrum 
suggest a density of 10 -40 Mpc -3 [3]. Therefore, we provide catalogs for the source 
density hypotheses 10" 35 Mpc" 3 , 10" 40 Mpc" 3 and 10" 45 Mpc" 3 . 

1. Mock Galaxy Catalogs 

To have a realistic treatment of CRs in the GZK transition region from 40-60 EeV, 
the galaxy catalog used to produce the UHECR mocks must be volume-limited to a 
much larger depth than is available in present-day "all sky" galaxy surveys, the best 
of which are 2MRS[3] and PSCz[S]. In particular, the galaxy catalog from which the 
source catalog is derived must be considerably denser than 10 -35 Mpc" 3 to simulate 
a scenario in which the source catalog has that density, and yet not every galaxy 
hosts a UHECR source. Therefore, we make use of the "LasDamas" mock galaxy 
catalogs [6] which are created using LCDM simulations (S1m = 0.25, £Ide = 0.75, 
H = 70 km/s/Mpc, eg = 0.8 and n s — 1) combined with an algorithm for populating 
the DM halos with galaxies, that is tuned to reproduce the clustering of galaxies as 
measured in the SDSS[7]. 
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Figure 1. Local density of galaxies in Carmen02_Mr21_Master (left) and 
Carmen27_Mr21_Master (right) in a slice of thickness 4000 km/s centered on the 
origin. 



We start with two LasDamas "Carmen" catalogs which are volume-limited for 
M r < —21 (includes all galaxies with an absolute r-band magnitude of -21 and 
brighter) to a redshift z — 0.165, and have a number density of 4 x 10 -4 Mpc~ 3 
or log(n) = —3.41 . They are: 

• Carmen02_Mr21_Master.dat - 522769 galaxies 

• Carmen27_Mr21_Master.dat - 520277 galaxies 

The local structure in the two catalogs can be seen in Fig. [T] The galaxies in 
each file are ordered by increasing redshift and each has a unique ID with the format 
mgC(catID)_(N) where catID is either 02 or 27. N is a unique integer that corresponds 
to the galaxy position in the redshift ordered list. 

We randomly subsample each catalog to obtain sets of catalogs at lower densities. 
We create one catalog with log(n)=-3.5 for each of Carmen02 and Carmcn27: 

• mockgalsCarmen02,3.5.dat 

• mockgalsCarmen27,3.5.dat 

We create 3 independent mock UHECR-source catalogs with log(n) = -4 for each 
of Carmen02 and Carmen27: 

• mockgalsCarmen02,4.0,R(i).dat where i=l,2,3 

• mockgalsCarmcn27,4.0,R(i).dat where i=l,2,3 

We create 3 independent mock UHECR-source catalogs with log(n) = -4.5 for 
each of Carmen02 and Carmen27: 

• mockgalsCarmen02,4.5,R(i).dat where i=l,2,3 

• mockgalsCarmen27,4.5,R(i).dat where i=l,2,3 
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Figure 2. The density of galaxies per redshift bin in each of the 6 Carmen 
catalogs with p = 10~ 4 - 5 Mpc~ 3 . 



In all of these source catalogs, the file format is RA Dec z id. The ID field gives 
the ID of the galaxy to allow comparison with the relevant master catalog. 

The density of galaxies per redshift bin is shown in Fig. [2] for redshifts within z = 
0.06. The error bars on the histogram show the Poisson error for the number of galaxies 
in the bin. The lack of any galaxies within z = 0.0045 for Carmen 02 and within z = 
0.0025 for Carmen 27 is a feature of the master catalogs. Carmen02_Mr21_Master.dat 
contains no galaxies within z = 0.0045, Carmen27_Mr21JVlaster.dat has 18 galaxies 
within z = 0.0045, but none within z = 0.0025. The lack of very nearby galaxies is 
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Figure 3. The number of galaxies within z ra dius of a candidate central galaxy, 
chosen at random from the master catalogs of Carmen 02 and Carmen 27. Each 
plot shows the results of 200 randomly chosen central galaxies. 



intentional, to simulate the underdense region around the Milky Way. 

Fig. [3] shows the result of taking 200 random galaxies from the full catalog 
(Carmen02_Mr21_Master.dat or Carmen27_Mr21_Master.dat) as central reference 
points, and counting the number of galaxies within a spherical volume with radius 
Zradius = 0.0025 or z ra di U s = 0.0045 around each central point. This shows that an 
underdense region within z ~ 0.0045 is not a generic feature for a randomly selected 
central galaxy. 

2. Mock Cosmic Rays 

From these mock-source catalogs, we create mock UHECR catalogs from a given mock- 
source catalog by choosing a galaxy at random, with a probability Pi — Wi/T,jWj, 




Figure 4. Distribution of energies at Earth, as a function of redshift of source. 



with iUj = rj 2 , where r» is the distance of the ith galaxy. (Distance, not redshift, 
is used here because the effects of peculiar motion can lead to up to 20 Mpc shifts 
in the distance with respect to that calculated from the Hubble flow.) A CR from 
this galaxy is then assigned an energy according to an injection energy spectrum and 
propagated to Earth. Energy loss from redshifting as well as photo-pion and e+e- 
production is included. In the simulations presented here, the CR is assumed to be a 
proton; the injection spectrum is proportional to E~ 2 Z and ranges from 40 to 1000 
EeV. If the energy after propagation is greater than 40 EeV, this proton is kept as 
an event. This procedure of source sampling and propagation is repeated until the 
desired number of events have been kept. Fig. [4] shows the final energy distribution as 
a function of redshift (actually, distance expressed as an equivalent redshift without 
peculiar motion). The resulting catalogs are: 

• mockgalsCarmen02,3.5,CRs.dat 

• mockgalsCarmcn27,3.5,CRs.dat 

• mockgalsCarmen02,4.0,R(i),CRs.dat, i=I,2,3 

• mockgalsCarmen27,4.0,R(i),CRs.dat, i=I,2,3 

• mockgalsCarmen02,4.5,R(i),CRs.dat, i=I,2,3 

• mockgalsCarmen27,4.5,R(i),CRs.dat, i=I,2,3 

Each cosmic ray catalog includes 100,000 mock cosmic rays, generated using the 
corresponding mock-galaxy catalog. The file format is: 

CRId E(EeV) srcID srcZ srcRA srcDec InjectionE RalDeg ... Ra8Deg Dec8Deg 

• CRId is an integer id 

• E is the energy of the CR at earth 

• srcID is the ID of the source galaxy 

• srcZ, srcRA and srcDec are the redshift, RA and Dec of the source galaxy 

• InjectionE is the energy of the cosmic ray at the source (in eV) 

• Ra(j)Deg, Dec(j)Deg give smeared arrival directions for a smearing of j degrees. 
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The columns headed Ra(j)Deg, Dec(j)Deg provide a range of arrival directions for 
the cosmic ray assuming a Gaussian deflection of the arrival direction: 

dP{6) _ 1 
dil ~ 27UT 2 V 2cr 

For this distribution a cone of radius 1.5a contains 68% of the smeared events. 
This angle is designated ogs- Ra(j)D e g, Dec(j)Deg in the header columns refer to a^s 
to label the arrival directions with increasing angular deflections. 

The full files can be downloaded from the project web page 



3. Attribution 



For proper attribution, we ask: 

(i) Any use must credit the mock UHECR catalog project and cite this paper and 
the LasDamas paper about the galaxy catalog from which it is derived [Hj. 

(ii) If used in a publication, inform us about your work before submission. 

(iii) Let us know that you downloaded the files, by sending an email to gf25@nyu.edu, 
so we can alert you if there are any problems discovered and about improvements 
implemented. 

By downloading the mocks, you implicitly agree to these terms. 



4. Future Work 



The Carmen mock catalogs will be extended in the near future to include mixed 
composition cosmic rays. The work will also be extended to provide mock cosmic ray 
catalogs for the known local distribution of galaxies from 2MRS. 
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